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Two low-dimensional magnetohydrodynamic models containing three velocity and three magnetic 
modes are described. One of them (nonhelical model) has zero kinetic and current helicity, while 
the other model (helical) has nonzero kinetic and current helicity. The velocity modes are forced in 
both these models. These low-dimensional models exhibit a dynamo transition at a critical forcing 
amplitude that depends on the Prandtl number. In the nonhelical model, dynamo exists only 
for magnetic Prandtl number beyond 1, while the helical model exhibits dynamo for all magnetic 
Prandtl number. Although the model is far from reproducing all the possible features of dynamo 
mechanisms, its simplicity allows a very detailed study and the observed dynamo transition is shown 
to bear similarities with recent numerical and experimental results. 
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I. INTRODUCTION 

The understanding of magnetic field generation, usu- 
ally referred to as the dynamo effect, in planets, stars, 
galaxies, and other astrophysical objects remains one of 
the major challenges in turbulence research. There are 
many observational results from the studies of the Sun, 
the Earth, and the galaxies [H [21 E]. Dynamo has also 
been observed recently in laboratory experiments [U |S] 
that have made the whole field very exciting. Numerical 
simulations [6j [7l [3 [9l [10] also give access to many useful 
insights into the physics of dynamo. However, the com- 
plete understanding of the dynamo mechanisms has not 
yet emerged. 

The two most important nondimensional parameters 
for the dynamo studies are the Reynolds number Re = 
UL/v and the magnetic Prandtl number Pm = ^/Vi 
where U and L are the large-scale velocity and the large 
length-scale of the system respectively, and v and rj are 
the kinematic viscosity and the magnetic diffusivity of 
the fiuid. Another nondimensional parameter used in 
this field is the magnetic Reynolds number Rcm, defined 
as UL/rj. Clearly Rem = Re Pnn hence only two among 
the above three parameters are independent. Note that 
galaxies, clusters, and the interstellar medium have large 
Pm, while stars, planets, and liquid sodium and mercury 
(fluids used in laboratory experiments) have small P„ 

mm- 

In a typical simulation, the magnetofluid is forced and 
the dynamo transition is considered to be observed when 
a nonzero magnetic field is sustained in the steady-state 
laminar solution or in the statistically stationary tur- 
bulent solution, depending on the regime. Typically, 
dynamos occur for forcing amplitudes beyond a critical 
value which also defines the critical Reynolds number Re'^ 
and the critical magnetic Reynolds number Re^ ■ One of 
the objectives of both the numerical simulations and the 
experiments [¥! 151 is the determination of this critical 



magnetic Reynolds number Re'^ . It has been found that 
Re^ depends on both the type of forcing and the Prandtl 
number (or Reynolds number) , yet the range of Re^ ob- 
served in the numerical simulations is from 10 to 500 for 
a wide range of P,„ (from 5 x lO"'^ to 2500). 

In recent magnetohydrodynamics (MHD) simulations, 
Schekochihin et al. [BJ |TT] applied nonhelical forcings 
and observed that the dynamo is active for a magnetic 
Prandtl number larger than a critical Prandtl number 
Pm that is around 1. For fluids with small Prandtl num- 
ber Pm {Pm < l)i numerical simulations [7, 8, 9J indicate 
that the dynamo can be produced using forcings having 
local helicity (the net helicity of the force could still be 
zero) . The range of the critical magnetic Reynolds num- 
ber in most of the simulations [71[Hl[n] is 10 to 500. Note 
that in the Von-Karman-Sodium (VKS) experiment, the 
critical magnetic Reynolds number is around 30. 

There are many attempts to understand the above ob- 
servations. For large Prandtl number, the resistive length 
scale is smaller than the viscous scale. For this regime, 
Schekochihin et al. [6', TP suggested that the growth rate 
of the magnetic field is higher in the small scales because 
stretching is faster at these scales. This kind of magnetic 
field excitation is referred to as small-scale turbulent dy- 
namo. For low Pm Stepanov and Plunian [12, argue for 
similar growth mechanism. Their results are based on 
shell model calculations. The numerical results of Iskakov 
et al. |7] however are not conclusive in this regard. The 
main arguments supporting these explanations are based 
on the inertial range (small-scales) properties of turbu- 
lence [HI HI]- In this paper, we present an alternate view- 
point. We show that a low-dimensional dynamical sys- 
tem containing only large scales properties of the fields, 
three velocity and three magnetic Fourier modes, is able 
to reproduce some of the above numerical results. For 
certain types of forcing of velocity, a dynamo transition 
is observed for Rem > R^m- These observations indi- 
cate that the large-scale eddies may also be responsible 
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for the dynamo excitation, and several important prop- 
erties can be derived from the dynamics of large-scale 
modes. These observations are consistent with earlier re- 
sults of MHD turbulence indicating that the large-scale 
velocity field provides a significant fraction of the energy 
(around 40%) contained in the large-scale magnetic field 
[T1[T1[T51[T5J[T71[TH]. 



II. DERIVATION OF LOW DIMENSIONAL 
MODELS 

The incompressible MHD equations are: 

dtu = -n(u, u) + n(b, b) + i^V^u + f - \/ptot, (1) 
atb = -n(u,b) + n(b,u) + KV2b, (2) 
V-u = V-b = 0, (3) 

where ptot is the sum of the hydrodynamic and magnetic 
pressures divided by the density. The bilinear operator 
is defined as n(ai, 32) = ai • Va2. Only periodic solution 
in a cubic box with linear dimension £ will be consid- 
ered. The smallest non-zero wave vector is thus given by 
ko = 2-K jl. The model is derived by projecting the MHD 
equations on the subspace S"^ spanned by a small num- 
ber of basis vectors that are compatible with the periodic 
boundary conditions and can be regarded as a subset of a 
complete basis. Only three vectors will be considered in 
this study, so that the projection of both the equations 
for u and b on these vectors leads to a six-dimensional 
system of equations. The three vectors are defined as 
follows: 



62 = 



\/2T7? 



\/2T7? 



— sin k^x cos k^z 
h sin kox sin kgz 

cosfco^^ sinfcoz 

ft, sin kou sin k^z 

— sinfco?/ cosfco^ 
cosfcoy sin/coz 

— sinfcoa; cosfco?/ cos2fcoz 

— coskox sinfcoy cos2fco2; 
cos k^x cos fcoy sin 2/co^^ 





(4) 



(5) 



—2 sin k^x cos k^y sin 2k^z 
sin k^x sin fcoj/ cos 2kf,z 



(6) 



They depend on a free parameter h that determines the 
amount of helicity that can be carried on by the basis 
vectors. For ft = 0, each mode is nonhclical. The three 
vectors are also divergence free V • = and are or- 
thonormal ■ e^) — Sa/j- The inner product of two 
arbitrary vectors vi and V2 is defined by: 

(vi -va) = ^ J^dxJ dyj dzvi{x,y,z) ■V2{x,y,z). (7) 

For ft = 0, the basis functions ei and 62 have field con- 
figurations in xz and yz respectively. 



The projection operator is noted P and its application 
on a vector v is defined by: 



P[v] = v< = t;i ei -I- V2 62 -t- V3 63 



(8) 



where Va — (v • e^). The projected part of v in 
is noted and the difference with the original vector 
is noted = v — v^. The projection of the velocity, 
magnetic, and force fields are expressed as follows. 

P[u] U< = (Ui ei +U2 62+ U3 63) u*, (9) 

F[b] = b< = (61 ei + &2 62 + 63 63) ^i^ (10) 
P[f] = f< = (/i ei + /2 62 + /3 63) (u*)^ ko. (11) 

where the Mq's, b^s and /q's are dimensionless real num- 
bers and u* = vko. The 63 component of the forcing is 
supposed to be zero. The system of dynamical equations 
for these variables can be derived from Eqs. (T][3l. For 
instance the projection of the equation for the velocity 
on the vector ei, {dtU ■ ei), yields the equation for u'l. 
The resulting system of equations is: 



Ul 

h 



r{h) (1 - ft2) (u2?/3 - ^2^3) - 2 wi + /i (12) 
r(ft) (1 + 3ft2) (U1U3 - fei&s) -2u2 + .f2 (13) 
-2r(ft,) (1 + ft2) {uiU2 - &162) - 6 ii3 + /3(14) 
r(ft) (1 + ft2) (U263 - h2U3) - 2 F„-i 61 (15) 
-r(ft) (1 + ft2) (U361 - 63U1) - 2 P;;! 62 (16) 
-2r(ft,) e {uib2 - 61U2) - 6 63- (17) 



where r(ft) = 2/ {{2 + h^) + 10ft2) . 

The structure of this model is, of course, reminiscent of 
the original MHD equations. In particular, it is interest- 
ing to study how the conservation of the ideal quadratic 
invariants of the Navier-Stokes and MHD equations can 
be expressed in terms of the variables {ua, The con- 
servation of the kinetic energy Ek by the Navier-Stokes 
equations is a direct consequence of 



(n(u, u) • u) = 



(18) 



Because this equality holds for any velocity, it is also true 
for and, in the absence of magnetic field, the nonlin- 
ear terms in the system of equations ( T2][T7 l conserve the 
kinetic energy P^ = {ul + u\+ '«3)/2 associated to u<. 
The cross helicity = (wibi + U262 + U3&3) and the 
total energy, the sum of the kinetic energy and the 
magnetic energy = {h\ + b\ + b\) /2, are conserved by 
the nonlinear terms for the same reason. However, the 
conservation of the kinetic helicity = (u • w) by the 
Navier-Stokes equations and the conservation the mag- 
netic helicity Hm = (b • a) by the MHD equations, where 
w = V X u is the vorticity and a is the vector potential 
(b — V X a) have no equivalent in the system ( 12p7 l. 
Indeed, the conservation of the kinetic helicity is a con- 
sequence of 



(n(u, u) 



UJ 







(19) 
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However, both the nonlinear term and the vorticity are 
not fully captured in 5^, even when computed from u^: 



solutions: 



P[n(u<,u<)] ^ n(u<,u<) 
P[V X u<] 7^ V X u< , 



(20) 
(21) 



As a consequence, the kinetic helicity carried on by 
is not conserved by the nonlinear term in the low- 
dimensional model and, for the same reason, the mag- 
netic helicity is not conserved either. For instance, 
the non-conservation of the kinetic helicity by the low- 
dimensional model can be expressed as, 



(n(u<,u<)<-(Vxu<)<) = -(n(u<,u<)>-(Vxu<)>) , 

(22) 

and the right hand side in this relation in not represented 
in the equations ( 12p7 1 . In the following section we will 
solve the truncated MHD equations ( 12][T7 l under spe- 
cial cases. 



III. ANALYSIS OF THE MODEL 



Despite its simplicity, the complete analysis of the sys- 
T2|l7| is quite difficult. 



tem of equations ( 
simple search for t 
algebraic equations 



Indeed, the 

le fixed points leads to very intricate 
We have thus focused our analysis 
on two limit cases corresponding to strictly nonhelical 
equations {h = 0) and the case h — I, both of which can 
be treated analytically. 



Nonhelical model 



The system of equations ( 12][l7 1 is thus first considered 
for h = 0. In that case the kinetic helicity captured in 
subspace 5^ vanishes ((u< • (V x u<)) = 0). In this first 
model, the forcing is chosen to be 



f=^(e, + e,) 



(23) 



Flmd ^± < 



Ml 



U2 



V2 



V2 



/± v//2-1152 
/T v//2-1152 



(24) 



L "3 



= -2^6 



Fluid B ■ 



Us 



(■■^^ - 12) 
s 

(s2 - 12)2 



(25) 



54 s2 



1 /3 

where s = (9//\/2 + 3^192-^9/2/2) . The re are two 
additional fixed points with nonzero magnetic fields: 



V2fP^ 
4(P™ + 1) 



U3 



2y/6 
P 







&1 = 









Fluid A 






\ MHD 


Fluid B ^ 



(26) 



Pm 



FIG. 1: Plot of critical force /cj as a function of Pm for the 
nonhelical model. Dynamo is excited for / > /^j, only for 
Pm > 1. 



where / = \fW^- Hence, /i = /a = //V2 and /s = 
0. The fixed points can then be computed analytically. 
Three fixed points correspond to a vanishing magnetic 
(61 = &2 = ^3 = 0) field and will be referred to as fiuid 



The stability of the above fixed points can be estab- 
lished by computing the eigenvalues of the stability ma- 
trix. After some tedious algebra, it can be shown that the 
above fixed points are stable in the three regions shown 
in Fig.nTa) and defined by the following simple equations 
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in the plane {Pm, f)- 



P — — 1 
/ = A, =24n/2 



Pm + 1 

p3/2 



(27) 
(28) 

(29) 



Since the velocity and magnetic field amplitudes Ua and 
ba must be real numbers, the solutions Fluid only 
exist for / > /^j. They are stable for P,^ < 1. The 
solution Fluid B is stable for Pm < 1 and / < /d, and 
for P,77 > 1 and / < /c2- It is unstable elsewhere. The 
solution MHD^ is stable for P„i > 1 and / > and 
is unstable elsewhere. In summary, we have only fluid 
solutions for Pm < 1, while an MHD solution is possi- 
ble only for P,7^ > 1 and / > fc2- The above solutions 
cover the entire (/, Pttj) parameter space. The system 
converges to one of these states depending on the pa- 
rameter values irrespective of its initial conditions, and 
the system is neither oscillatory nor chaotic. Clearly, the 
ampUtude of bi and 62 close to the dynamo threshold 
increases as \/ f — fc2- The dynamo transition in this 
low-dimensional model is thus a pitchfork bifurcation. 

The above results compare remarkably well with the 
numerical findings of Schekochihin et al. [SJ [TT] where a 
dynamo transition is also found for P,7i > P^ with P^ 
near 1. Considering the simplicity of the above model, 
such a quantitative agreement may very well be fortu- 
itous. Nevertheless, it is also interesting to notice that, 
in our model, /cj decreases for increasing Prandtl num- 
bers. Hence, it is easier to excite nonhelical dynamo for 
larger P„j or more conductive magnetofluid. This feature 
of our model is also in agreement with recent numerical 
simulations [6l [7l 13 [9l [10] where an increase of the critical 
magnetic Reynolds number Pe^ is reported for smaller 
values of Pm ■ 

It is also interesting to express the above results in 
terms of the Reynolds number instead of the forcing am- 
plitude. A Reynolds number can be build by defining the 
large scale velocity as Ul = \/u'{ + u\ u*. Based on this 
velocity scale, the kinetic Reynolds number is given by: 



Re^ — = < ^ 



f Pn 



I 2(1 + P„7) 



Fluid A 



Fluid B 



MHD 



(30) 



For large /, the amplitude of the magnetic field is pro- 
portional to the kinetic Reynolds number: 



Re I 



(31) 



This result is consistent with the result obtained by Mon- 
chaux et al. [ll[51[Tn]. The critical Reynolds number for 



the dynamo transition between the Fluid B and MHD 
solutions is easily computed: 



Re" = 6a 



and the critical magnetic Reynolds number Re'^ 



Rel^ ^PmRe^ Q^2P~„. 



(32) 



(33) 



Since P777 > 1, Pe^„ > 6%/2. Note that Re^Re" 72. 
Hence, the Re'^ — Rem curve is a hyperbola for Re'^ < 6 a/2 
since dynamo exists for P777 > 1. 

To gain further insights into the dynamo mechanism 
we have investigated the energy exchanges among the 
modes of the model. The fluid modes ui and U2 gain en- 
ergy from the forcing f and give energy to the mode U3. 
The magnetic modes hi and 62 gain energy from the mode 
M3. The net energy transfer to the mode 63 through non- 
linear interaction vanishes. Consequently the mode 63 
goes to zero due to dissipation. The 61 and 62 modes also 
lose energy due to dissipation. The energy balances for 
these modes reveal interesting feature of dynamo transi- 
tion. At the steady-state, the energy input in the mode hi 
due to nonlinearity {—bib2U3/V6 — —b^u s/VG) matches 
the Joule dissipation {2bl/Pm) [Eq. For f < f^, 

the value of U3 is less than —2\/6/Pm, hence bi ^ 
asymptotically, thus shutting down the dynamo. For 
/ > f^^, U3 = — 2^/6/ Pm, thus the energy input to the 
mode 61 exactly matches with the Joule dissipation. This 
is the reason why bi and 62 are constants asymptotically 
in the MHD regime. 

We also studied a variant of the above model in which 
fi — and /2 = /. For this model, the only solution is 
U2 — f /2 with all the other variables being zero. Hence 
this model does not exhibit dynamo. In the next subsec- 
tion we will discuss another low-dimensional model that 
contains helicity. 



B. Helical model 

In this second version of the model, the value h = 1 
has been chosen. The helicity captured in the subspace 
5^ is then given by: 

(u<.(Vxu<))^ Q {uj-uD-^ul^ {u*fk,. (34) 

If we restrict again the forcing to the large-scale modes 
(/3 = 0), the amount of helicity carried on by the velocity 
fleld is expected to increase, at least in absolute value, if 
the forcing act differently on the modes ei and 62. We 
have chosen the extreme case for which the forcing acts 
only on U2 (/i = /3 = and /2 = /): 



f =/e2 



(35) 
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where again / = •\/(f • f). For these parameters, a unique 
fluid stationary solution and two stationary MHD solu- 
tions are found: 



ui = M3 = 



Fluid ■ 



(36) 



' ui = U3 = 



U2 



6 V3 
p 



72 
p 



2V3f 



1/2 



(37) 

All these solution obviously carry a non-zero resolved he- 
licity as defined by (34 1. The fluid solution is stable for 
/ < /c3 = 12 a/3 /Pm, while the MHD solutions are sta- 
ble for f > fc3- The plot of the critical force fc^ vs. Pm 
for helical model is shown in Fig. |2] Note that the helical 
model exhibits dynamo for all P„i as long as f > fc^. 
Also, the critical forcing for the helical model is lower 
than the corresponding value for the nonhelical model. 
Hence it is easier to excite helical dynamo compared to 
the nonhelical dynamo consistent with recent numerical 
simulations [3 [i [i [TD] . 




Pm 



FIG. 2: Plot of critical force amplitude fc^ as a function of 
Pm for the helical model. The dashed lines reproduce the 
stability regions of the non-helical model. 



Here again, the Reynolds number based on the large 
scale velocity Ul = u* y/u1 + U2 can be defined and 



yields: 



Re = 



lyko 



f/2 Fluid 
6\/3 



(38) 



P„ 



MHD 



and the magnetic Reynolds number for MHD case is 
R^m = 6-\/3, a constant. 



Using Eq. (34 1 we can compute the resolved kinetic 

ifined a 



helicity Hk and current helicity Hj defined as 
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(u< • (V X u<)) 
144, 



H.J = 




iu*f ko 



(39) 



(40) 



Pouquet et al. [20], and Chou |2T] conjectured the alpha 
parameter of dynamo to be of the form 



1 



at 



r (-U • V X u-fb • V X b) , (41) 



where r is the velocity de-correlation time. Clearly a 
is optimal if Hk < and Hj > 0. Similar features 
are observed in fiux calculation of Verma [22]. These 
conditions are satisfied in the helical model indicating a 
certain internal consistency with the results of Pouquet 
et al. [20], and Chou [21]. 

The energy exchange calculation of the helical model 
reveals that the magnetic modes 61 and 63 receive energy 
from the velocity mode U2; the rate of energy transfer 
is proportional to &1&3M2 that matches exactly with the 
Joule dissipation rate. Since U2 — &^/?,/Pm for all 
the rate of energy transfer 6163U2 matches with Joule dis- 
sipation rate (oc h\/Pjn), and the dynamo is possible for 
all Pm in case of helical model. In contrast, in nonheli- 
cal model the corresponding velocity mode U3 varies as 
for Pm > 1, but saturates at 2\/3 for Pm < 1, 
hence the energy transfer cannot match the Joule dissi- 
pation rate for Pm < 1 for nonhelical model thus shutting 
off the dynamo for Pm < 1. This is one of the main dif- 
ference in helical and nonhelical models. Also note that 
as discussed at the end of the previous subsection, the 
nonhelical model with the forcing /i = 0, /2 = / does 
not exhibit dynamo for any parameter. 



IV. DISCUSSION 

In this paper, a class of low-dimensional models that 
exhibit dynamo transition is derived. These models de- 
pend on several parameters such as the Prandtl number, 
the forcing amplitude and a parameter, h, that charac- 
terizes the ability of the velocity modes to carry kinetic 
helicity. The fixed points of two simple models corre- 
sponding to ft, = and /i = 1 are studied in details. The 
first model is nonhelical (zero kinetic and current helic- 
ity) and is compatible with a stationary nonzero mag- 
netic solution only for Pm > 1. The second model has 
nonzero kinetic and current helicities, and it has nonzero 
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stationary magnetic solution for all Pm- These findings 
confirm the idea that both kinetic and current helicities 
may play an important role in the dynamo transition, 
especially in helical model for Pm < 1- These two val- 
ues of h correspond to the only values for which one of 
the nonlinear coupling in the low dimensional model van- 
ishes. With the specific choice of the forcing proposed in 
the previous section, these models have then exact and 
tractable solutions. Arbitrary values of h would lead to 
much more complex systems. 

Obviously, the models presented here are only a small 
subset of many possible MHD models that could exhibit 
dynamo. For instance, Rikitake |,23j constructed a dy- 
namo model consisting of two-coupled disks that shows 
field reversal. Nozieres [57 proposed a model involv- 
ing one velocity and two magnetic field variables. These 
two models are phenomenological. Recently, Donner et 
al. [25] proposed a truncation of the MHD equations 
along the same lines as our model, but derived a much 
more complex system of equations for 152 modes. Don- 
ner et al. |25j analyzed the dynamical evolution of these 
modes for P^ = 1 and observed steady-state and chaos 
in their system. The main advantage of our model is that 
it allows a complete analytical treatment. 

The six variables of our model are only representative 
of the large-scale modes of the system, while a realis- 
tic description of a turbulent system exhibiting dynamo 
should have a large number of modes. Although the 
small scale variables are often quite important in tur- 
bulent flow, in some cases the large scale variables may 
determine some of the flow properties. The surprising 
success of the very simple models proposed here in re- 
producing several feature of the dynamo transition could 
suggest that we are in such a situation. Schekochihin 
et al. [5J[TT], Stepanov and Plunian [T^, and Iskakov et 
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